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By employing the first-principles GW-Bethe-Salpeter Equation (BSE) simulation, we obtain, for 
the first time, the accurate quasiparticle (QP) band gap, optical absorption spectra and their de- 
pendence on the gate field of gated bilayer graphene (GBLG). Many-electron effects are shown to be 
extremely important to understand these excited-state properties; enhanced electron-electron inter- 
actions dramatically enlarge the QP band gap; infrared optical absorption spectra are dictated by 
bright bound excitons. Our results well explain recent experiments and satisfyingly solve the puzzle 
about the inconsistency between experimentally measured transport and optical band gaps. More- 
over, our calculation reveals fine excitonic structures and predicts exotic excitonic effects that have 
not been observed yet, which can be of interest for optoelectronics applications based on GBLG. 

PACS numbers: 73.22.Pr, 78.67.Wj, 71.35.Cc 



Despite its intriguing electronic, thermal and op- 
tical properties [l|-[3j, intrinsic graphene is a gapless 
semimetal, impeding its utility in bipolar devices and 
subsequent broad applications. Therefore, huge efforts 
have been made to overcome this barrier by generating a 
finite band gap in graphene or its derivatives 0-0- One 
promising approach is to apply the gate electric field per- 
pendicular to the AB (Bernal) stacked bilayer graphene 
?BLG) to break the inversion symmetry of sublattices 
0, [H| . Such an induced band gap of GBLG can be tuned 
in a wide range by field strength [H, [lll - [l3| , offering an 
important degree of freedom to optimize performance of 
graphene devices. 

However, we are still in lack of a satisfying under- 
standing on fundamental properties of GBLG, such as 
its quasiparticle (QP) band gap and optical excitations. 
For instance, electrical conductance experiments pj]-[l3j 
have confirmed the existence of a finite QP band gap but 
their measured value is annoyed by the inevitable contact 
resistance between electrodes and graphene sheet; while 
noncontacting optical measurements j&LL0| have revealed 
a tunable band gap in GBLG but these results are in- 
direct because the optical absorption peak (edge) is not 
conceptually equivalent to the QP band gap [LjUil- Par- 
ticularly, enhanced excitonic effects often dramatically 
shift the optical absorption peak as we have seen in many 
other reduced dimensional semiconductors [H, Qj) , mak- 
ing this inconsistence even more serious. Therefore, an 
accurate calculation with many-electron effects is urgent 
to settle the above inconsistency. 

Conventional density functional theory (DFT) simu- 
lations cannot answer above questions because of their 
known deficiencies of handling excited-state properties 
BES Est Tight-binding models (l§ have revealed ap- 
pealing physics of excitons in GBLG, but it has to rely 
on parameters. In particular, recent ab initio GW-BSE 
simulation has successfully predicted enhanced many- 
electron effects on intrinsic graphen e I20|, which are con- 
firmed by subsequent experiments [2lM23j. Therefore, a 
reliable first-principles calculation with many-electron ef- 
fects included is also promising. 



More importantly, beyond providing reliable parame- 
ters for device design, obtaining knowledge about excited 
states of GBLG will be of fundamental interest because 
it will fill our knowledge gap on many-electron interac- 
tions in two-dimension (2-D) narrow-gap semiconductors, 
a field that has not been well understood yet. In fact, 
it is challenging for first-principles simulations to accu- 
rately capture the nearly metallic electronic screening of 
narrow-gap semiconductors. For this purpose, an im- 
proved algorithm has to be developed and shall be of 
broad interest for the electronic-structure community. 

In this Letter, with the modified model accurately de- 
scribing the screening, we conclude three important re- 
marks about excited states of GBLG. 1) The QP band 
gap and its dependence on the gate field are obtained. 
The self-energy correction is significant because of the 
enhanced electron-electron (e-e) interactions; the calcu- 
lated QP band gajDS are more than 150 % of previous 
DFT predictions [24], which is beneficial for device 
applications since a wider band gap implies higher work- 
ing temperature. 2) Optical absorption spectra of GBLG 
are dominated by excitonic effects. With electron- hole (e- 
h) interactions included, our calculated absorption peaks 
are in excellent agreement with recent experiments [ij], 
well explaining the inconsistency between QP gap and 
optical gap. 3) Excitons in GBLG exhibit a number of 
unusual features. For example, the electron and hole 
of the lowest-energy dark exciton are completely sepa- 
rated onto two graphene layers, giving rise to an optical 
approach to polarize BLG. Moreover, this separation of 
electron and hole offers a neat opportunity to evaluate 
entangling effects, such as the exchange interaction, of 
many-electron systems. 

To reveal the significance of many-body correlations 
in GBLG, we perform the calculation of excited states 
as the following procedure. First, ground-state energy 
and wavefunctions are obtained by DFT within the local 
density approximation (LDA). Secondly, the QP energy 
is calculated within the single-shot GW approximation 
[3. Finally we obtain the exciton energy, wavefunctions 
and optical absorption spectra by solving the following 
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where A^ c ^ is the amplitude of excitonic state S consist- 
ing of single-particle hole state \v\s) and electron state 
\ck). K eh is the e-h interaction kernel and ft s is the ex- 
citon excitation energy. E c k and E v k are QP energy of 
electrons and holes, respectively. 

All calculations are based on plane-wave basis and 
norm-conserving pseudopotential approximations with a 
60-Ry energy cutoff. To eliminate the spurious inter- 
action between neighboring BLG, the slab-truncation 
scheme is applied to mimic isolated GBLG [26|, |27j. The 
electric field is applied by the periodic sawtooth potential 
perpendicular to graphene layers. 

The crucial part of describing many-electron interac- 
tions is to obtain the dielectric function. For GBLG with 
the truncated Coulomb interaction, the inverse static di- 
electric function e _1 (q) rapidly changes within the long 
wave-length regime q — >• 0, which is similar to what have 
been noticed in recent first-principles simulations of car- 
bon nanotubes (CNTs) [l6j|. A brute- force way to cap- 
ture this feature is to use a dense q-grid, which demands 
formidable computational resources. To solve this prob- 
lem, we deliberately employ the mini Brillouin zone (BZ) 
sampling scheme to count in this sharply-varying char- 
acter as motivated by Refs. [26|, l27j, and extend it to 
both evaluating the QP energies and solving the BSE 
(see the supplementary material for details). As a re- 
sult, a 72 x 72 x 1 coarse k-grid sampling is adequate for 
the GW calculation. In addition, we employ a partial 
1440 x 1440 x 1 fine k-grid sampling around the Dirac 
cone for a dependable BSE calculation. 

The LDA and GW band structures near the BZ cor- 
ner (the K point) are plotted in Fig. [1] (a) for GBLG, 
respectively. The applied gate field induces a finite band 
gap and changes the band dispersion to the Mexican-hat 
feature. After including the self-energy correction via 
the GW calculation, the Mexican-hat-shaped feature re- 
mains intact, nevertheless the fundamental band gap is 
significantly enlarged due to the depressed screening of 
isolated GBLG. Moreover, the slope of band dispersion is 
sharpen by the self-energy correction, implying a smaller 
effective mass of free carriers. 

We also investigate the QP band gap dependence on 
the applied gate field as shown in Fig. [1] (b). The QP 
band gap can be varied from zero up to 300 meV under 
experimentally reachable gate field, which is also more 
than 150% of previous DFT predictions. This larger tun- 
able range and a wider gap are desired for device appli- 
cations because a wider gap means higher working tem- 
perature and lower noise. Moreover, if listing the ratio 
of the self-energy correction to their DFT/LDA value, 
we see it is 56%, 67%, 78% and 81%, respectively, as 
the applied field is decreased. This growing trend of the 
correction ratio for the smaller gap is of particular in- 
terest because recent experiments [13] shows a possible 
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FIG. 1: (Color online) (a) DFT/LDA and GW calculated 
band structures around the Dirac point of BLG under a gate 
field of 2 V/nm. (b) Comparison of the "gap" values ob- 
tained from different approaches and their dependence on the 
gate electric field. The value of the optical gap is defined by 
the position of the first bright peak of the optical absorption 
spectrum. The experimental values are extracted from Ref.0 



small band gap (around a few meV) even for BLG in the 
absence of gate field. However, due to the limit of the 
accuracy of our numerical simulation, we cannot resolve 
those energy differences below 10 meV and hence more 
advanced simulation technique needs to be developed. 

In addition, the recent optical measurements of the op- 
tical gap are plotted in Fig. ffl(b) as well. The key feature 
is that the QP gap is substantially larger than both pre- 
vious DFT predictions I24L [25| and measurements from 
optical experiments, [9|, |28|. The inconsistency QP band 
gap with optical measurements has also been observed 
in several other semiconducting nanostructures [H, [I?} , 
which manifests enhanced excitonic effects and motivates 
the following calculation on optical spectra of GBLG. 

Fig. [2] depicts the optical absorption spectrum and its 
evolution subject to the increasing field magnitude. We 
first focus on absorption spectra in the absence of the 
e-h interactions (blue lines in Figs. 12(a)). In the low- 
energy regime, the absorption is mostly contributed by 
the transition from the highest valence band (vl) and the 
lowest conduction band (cl). As expected, the absorp- 
tion onset displays a blueshift as the electric field uplifts 
the band gap magnitude. Meanwhile, the prominent ab- 
sorption feature is gradually broadened and split into a 
double-peak structure [I\ and I2) which stems from the 
two one-dimensional-like von-Hove singularities [HI, [29| 
at opposite "Mexican-hat brims" (Fig. [T] (a)), which is 
consistent with previous DFT results [25|. 

Surprisingly, the von-Hove singularity at the K point, 
however, does not contribute greatly to the absorption 
and therefore is not resolved in the spectra. This is be- 
cause the relevant valence state |vk) and the conduction 
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FIG. 2: (Color online) (a) Optical absorption spectra of 
GBLG. The vertical black dashed line indicates the GW fun- 
damental gap. The incident light is polarized parallel to the 
graphene plane. A 10 meV Gaussian smearing is applied, (b) 
Optical activity and eigenenergy of excitons. Each bar rep- 
resents one exciton state and its height (plotted in the log- 
arithmic scale) indicates the corresponding optical activity. 
The lowest-energy dark exciton D and the prominent exciton 
A are particularly outlined by widened dark and red bars, 
respectively. 



state |ck) at the Dirac point K are strongly localized on 
different layers upon field-induced symmetry breaking. 
Therefore, the overlap of wavefunctions is very small, 
which eventually leads to a negligible oscillator strength. 
This can be seen in Fig. [3] (a), in which we present the 
contour plot of the oscillator strength around the cor- 
ner of the first BZ. The strongest oscillator strength is 
actually from the "Mexican- hat brims" regime while it is 
almost zero at the K point. 

With e-h interactions included, a different optical ab- 
sorption picture emerges. As shown in Figs. [2] (a), the ex- 
citon effect dramatically reshapes the spectra; the broad 
asymmetric absorption peak in the single-particle picture 
is replaced by a symmetric prominent absorption peak. 
This peak lies below the QP band gap, indicating the 
existence of bound e-h pairs. The binding energy under 
four sampling voltages is 35, 54, 76 and 80 meV, respec- 
tively, which is also close to previous tight-binding cal- 
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FIG. 3: (Color online) (a) The distribution of single-particle 
oscillator strength in the reciprocal space. We only include 
transitions from the highest valence band to the lowest con- 
duction band, (b) and (c) The distributions of the square of 
the exciton amplitude (|A^ ck | 2 ) of the dark exciton D and 
bright exciton A in the reciprocal space. The square black 
dots mark the three identical locations of the minimum en- 
ergy gap. 



culations [19]. Remarkably these peak-position is in ex- 
cellent agreement with the previous infrared microspec- 
troscopy experiment [9| as shown in Fig. [1] (b). Under 
realistic experimental setups, both self-energy corrections 
and e-h interactions shall be reduced by the screening ef- 
fect of dielectric substrates. On the other hand, these 
reductions may cancel each other more or less [30|, re- 
sulting in such a good match of our calculations with 
experimental data. 

In the higher energy regime (around 0.4 eV) next to 
the first optical active peak, the absorbance maintains 
a constant on the whole (~ 3%), which is significantly 
smaller than 4.6%, the ideal value of the optical ab- 
sorbance of BLG [32|, [331 ] - This is due to the sum rule 
of oscillator strength [14[ because e-h interactions drains 
the absorbance from the high-energy regime to enhance 
the exciton peak. 

It has to be pointed out that electron-phonon coupling 
shall be another important factor in determining the in- 
frared optical spectra of GBLG. For example, a G-mode 
phonon at 195meV has been found to be in Fano interfer- 
ence coupled with e-h excitations in GBLG [3l|. There- 
fore, we may expect this dip feature from such a G-mode 
phonon may impact the lineshape of our studied exciton 
absorption peaks. 

A close inspection of solutions of the BSE reveals an 
intriguing exciton picture that have not been observed by 
experiments. We plot the oscillator strength of excitons 
in a logarithmic scale in Figs. [2](b). The isolated exciton 
state with the largest oscillator strength, A, is responsi- 
ble for the symmetric prominent absorption peaks in the 
spectra. Surprisingly, there is one lower-lying excion, D, 
with a much weaker oscillator strength for most gated 
fields (except 4V/nm). This is contrary to the usual 
effective-mass model, in which the lowest singlet exciton 
shall be the brightest one involved with two bands. 

Furthermore, we observe that both the position and 
oscillator strength of this dark exciton D are more sensi- 
tive to the gate field than those of bright exciton A. As 
plotted in Figs. [21(b), the energy of D progressively ap- 
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proaches that of A with an increasing gate field strength 
and its optical activity is strongly quenched simultane- 
ously. In particular, when the gate field is more than 
3V/nm, the order of the bright and dark excitons is 
switched as shown in Fig. [T] (b). This tunable energy 
difference can surely affect the thermal population of ex- 
citon states and their luminescent performance. Such a 
tunable energy order of excitons is in qualitative agree- 
ment with previous tight-binding studies [19| . 

To understand the brightness of these exciton states, 
we need to further investigate the origin of their optical 
activity. For a typical field strength (2 V/nm), Figs. 
[3] (b) and (c) display the distribution of the square of 
exciton amplitude A^ ck for those two interested excitons 
A and D. Since the optical activity of an exciton i [l4| 
is 

\(0\v\i)\ 2 = \J2Kck(vmck)\ 2 , (2) 

vck 

which roughly means the product of the single-particle 
oscillator strength shown in Fig. [3] (a) and exciton am- 
plitude shown in Figs. [3] (b) or (c), we immediately see 
the product of exciton D is much bigger than exciton A, 
suggesting their markedly different brightness. 

Fig. 2] visualizes exciton wavefunctions in the real 
space. As is readily seen, both excitons A and D are 
strong charge transfer excitons but with distinct charac- 
ters. In particular, the electron and hole of the dark ex- 
citon D almost become disentangled. As shown in Figs. 
HI (a), (c) and (d), the electron and hole wavefunctions of 
exciton D are nearly completely separated into two layer. 
This is very dissimilar to the e-h correlation in other 2-D 
semiconductors [2O, [Hf . From the perspective of opto- 
electronic application, exciton D could yield interesting 
possibility of efficient e-h separation and polarize BLG 
by optical excitations. For the exciton A, the degree of 
e-h separation is much lower. In Figs. H] (b), (e) and 
(f), the electron distributes over a ring on the top layer 
while on the bottom layer the electron distributes on a 
disk centered at the hole. 

Moreover, these excitonic wavefunctions will be cru- 
cial to understand why the dark exciton D and the dark 
exciton A respond very differently to the electric field. 
As concluded in Fig. [1] (b), the energy level of exciton 
D exhibits an approximately linear relationship with the 
field strength, whereas that of exciton A shows a nonlin- 
ear behavior. This can be rationalized by the fact that 
exciton D can be viewed as a plane of dipole composed 
of dissociated electron and hole as revealed in Fig HI (a), 
whose energy levels of the positive and negative poles lin- 
early depend on the applied gate field. In contrast, the 
electron and hole for the exciton A are spatially entan- 
gled and therefore the energy level is less sensitive to the 
gate field and does not follow a simple linear trend. This 
explains the origin of the energy order switch when the 



applied gate field is more than 3 V/nm. 

In conclusion, we have provided first-principles calcu- 
lations for the QP energy and excitonic effects of GBLG. 
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FIG. 4: (Color online) (a) and (b) Side views of the isosurface 
plot of the square of wavefunctions of the excitons D and A. 
(c) to (f) Top view of these exciton wavefunctions for top and 
bottom layers, respectively. The hole position is marked by 
the open circle in (a) and (b) while it is located at the center 
of the bottom layer in (d) and (f). 



E-e and e-h interactions are significant and must be con- 
sidered to understand the electronic structure and opti- 
cal excitations of GBLG. Moreover, our calculation well 
explains recent experiments and reveals more physics as- 
sociated with many-electron effects. Finally, we have ob- 
served an exotic dark exciton structure that is not likely 
to present in conventional direct band gap semiconduc- 
tors. The different degree of charge transfer for different 
exciton states may be useful in optoelectronic applica- 
tions. 

The computational resources have been provided by 
Lonestar of Teragrid at the Texas Advanced Comput- 
ing Center. The ground state calculation is performed 
by the Quantum Espresso. The GW-BSE calculation is 
done with the BerkeleyGW package [35| with minor mod- 
ifications for excited states calculations in GBLG. 
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